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Unruh-DeWitt detectors interacting locally with a quantum field are systems under consideration 
for relativistic quantum information processing. In most works, the detectors are assumed to be 
point-like and therefore, couple with the same strength to all modes of the field spectrum. We 
propose the use of a more realistic detector model where the detector has a finite size conveniently 
tailored by a spatial profile. We design a spatial profile such that the detector, when inertial, 
naturally couples to a peaked distribution of Minkowski modes. In the uniformly accelerated case, 
the detector couples to a peaked distribution of Rindler modes. Such distributions are of special 
interest in the analysis of entanglement in non-inetial frames. We use our detector model to show 
that the noise detected in the Minkowski vacuum and in single particle states is a function of the 
detector's acceleration. 
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I. INTRODUCTION 

Relativistic quantum information is a multifaceted 
area that incorporates key principles from quantum field 
theory, quantum information and quantum optics to an- 
swer, primarily, questions of an information theoretic na- 
ture. In order to implement quantum information tasks 
in relativistic settings it is necessary to find suitable lo- 
calised systems to store information. Moving point-like 
detectors coupled to quantum fields have been consid- 
ered to carry quantum information in spacetime [IH3], 
perform teleportation [3] , and extract entanglement from 
the Minkowski vacuum [2j [5j |6] . Our research program 
aims at developing new detector models which are more 
realistic and simpler to treat mathematically so that they 
can be used in relativistic quantum information process- 
ing. 

In this paper we propose the use of finite-size detectors 
[THS] which are not only more realistic but also have the 
advantage of coupling to peaked distributions of modes. 
We design gaussian-type spatial profiles such that an uni- 
formly accelerated detector naturally couples to peaked 
distributions of Rindler modes. By expanding the field in 
term of Unruh modes, we show that the accelerated de- 
tector couples simultaneously to two peaked distribution 
of modes corresponding to left and right Unruh modes. 
As expected, the same detector interacts with a gaussian 
distribution of Minkowski modes when it follows an iner- 
tial trajectory. In the Minkowski vacuum, the response 
of the detector has a thermal signature when it is uni- 
formly accelerated and the temperature depends on the 
proper acceleration of the detector. 

In the prototypical studies of quantum entanglement 
in non-inertial frames, observers are assumed to analyse 
states involving sharp frequency modes [TOl [H] . In par- 
ticular, recent works analysing the entanglement degra- 
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dation between global modes seen by uniformly acceler- 
ated observers consider states of modes labelled by Un- 
ruh frequencies |12H14j . Our analysis provides further 
insight into the physical interpretation of the particle 
states which were analyzed in these works. Finite-size 
detectors are suitable to discuss such results from an op- 
erational perspective since sharp frequency modes are an 
idealization of the peaked distributions finite-size detec- 
tors couple to. We show that the response of the finite- 
size detector when the state of the field has a single par- 
ticle labelled by a Unruh frequency has a thermal term 
plus a correction with noise that depends on acceleration. 
Therefore, a degradation of global mode entanglement in 
non-inertial frames as a function of acceleration should 
be detected by uniformly accelerated finite-size detectors. 
Although global mode entanglement cannot be detected 
directly it can be extracted by Unruh-DeWitt type de- 
tectors becoming useful for quantum information tasks. 
Throughout our work we use natural units and the metric 
has signature (-, +, +, +). 

The action of a detector interacting with a quantum 
field in the interaction picture is given by [7J |S] , 

Si = J dVM-fi, (1) 

where dV is the volume element for the spacetime and M 
is the monopole moment of the detector parametrised by 
arbitrary coordinates in a (3 + l)-dimensional spacetime. 
We will assume a flat spacetime and the field <f> to be a 
real massive scalar field which satisfies the Klein-Gordon 
equation, 

V M V M - m 2 4> = 0, (2) 

where m > is the free field mass and V M denotes the 
covariant derivative on the spacetime. We require the 
detector to have, in a comoving reference frame (r, £), 
a constant spatial profile. In this case the monopole 
moment of the detector factorizes into a temporal op- 
erator valued function and a spatial function M(t, £) = 
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M(t)/(C) [T5lU7| . The temporal operator valued func- 
tion describes the internal structure of the detector which 
we model by a two-level system of characteristic energy 
A therefore, M(r) := cr_e~ iAT + a\e lAr . The operators 
cr + and <7_ create and annihilate, respectively, excitations 
in the internal structure of the detector. In this case, the 
action can be written in terms of the detector's proper 
time r, 

Si = f cItM(t)4>(t), (3) 

where the field the detector couples to is given by, 

0(r):= f d 3 (^J)f(()cb(T,0 (4) 

and det((; AI j y ) = g is the determinant of the metric tensor. 
We assume that the center of the detector follows a clas- 
sical trajectory in spacetime and the spatial profile /(C) 
determines how the detector couples to the field along the 
trajectory. This function, which must be real to ensure 
that the action is Hermitian, can be interpreted as a po- 
sition dependent coupling strength. Consider uc,(£(t)) 
to be field solutions to the Klein-Gordon equation evalu- 
ated along a point-like worldlinc parametrised by r cor- 
responding to the center point of the detector. The fre- 
quencies v of the modes are determined by an observer 
comoving with the center of the detector. The Hamilto- 
nian in terms of these modes takes the form, 

Hi = M(t)- J d 3 i)f(D)(u i> (T)a e + h.c.), (5) 

where at and ap are creation and annihilation opera- 
tors associated to the field modes of frequency v. The 
frequency distribution f(y) corresponds to a transfor- 
mation of /(C) into frequency space. In the ideal case 
where the detector is considered to be point-like, the 
spatial profile is /(C) = <S 3 (C " CO (here 8 3 (( - (') ■= 
— Ci)^(C2 — C^^CCa - C3) is the three dimensional Dirac 
delta distribution). The detector couples locally to the 
field and the coupling strength is uniformly equal for all 
frequency modes. When we model a finite-size detector, 
which corresponds to a more realistic situation, the de- 
tector couples naturally to a distribution of field modes. 
The frequency distribution will be determined by the spa- 
tial profile. In this sense the field 4>(t) corresponds to a 
window of frequencies. 

II. INERTIAL TRAJECTORY 

Now we specify a trajectory for the detector. When 
the center of the detector follows an inertial trajectory it 
is convenient to use Minkowski coordinates (t,x) where 
x := (x, y, z). In this case, the proper time of a comoving 
observer is r = t and we can also write the comoving 
spatial coordinates as C, = i. The solutions to the Klein- 
Gordon equation correspond to plane waves, 



where the frequency of the mode uj = Vfc ■ k + m 2 is 
strictly positive and k e M 3 denotes the momentum 
k {k x ,ky,k z ) . The inner product satisfies (uj:,up) = 
5 3 (k - k') [IS]. In this case the creation and annihilation 
operators associated to the Minkowski field modes sat- 
isfy [ a fc: a jj,,] = S 3 (k - k'). The field can be expanded in 
Minkowski modes as 

(j){t,x) = J d 3 k(a j :u i: (t,x) +h.c.) (7) 

From this, the frequency distribution expressed in 
Minkowski modes is, 

f(k) = f <r.rf(.r)<- dl . (8) 

which is the Fourier transform of the spatial profile f(x). 

We now design a spatial profile tailored so that the 
corresponding frequency detection window of the detec- 
tor is a gaussian distribution of modes peaked around a 
Minkowski frequency A. This choice is motivated by early 
works on relativistic entanglement where the states anal- 
ysed involved sharp frequencies f2 and Q' . For example, 
the Bell state, 

|0) = -^=(|O)o|O)o' + |l)o|l)oO, (9) 

which was analysed in a flat (l+l)-dimensional space in 
the charged [T3] and uncharged massless bosonic case 
|12) . Entanglement for Bell states in non- inertial frames 
was also discussed for Dirac fields [Ml H9l [20] . Sharp 
frequency states, |l)n = ^olO), are an idealisation of our 
gaussian wave-packets of the form 

|1 A > = J dk$(\,k)a\\Q), (10) 

where $(A,/j) is a gaussian distribution centred around 
A [12 ■ Our detector model will be useful to investigate 
questions of entanglement in non-inertial frames from an 
operational perspective and extract entanglement for rel- 
ativistic quantum information processing. An interesting 
question, which we intend to address in following work, 
is how much entanglement can be extracted by our de- 
tectors from field mode entangled states as a function of 
the detector's acceleration. 

We find that a gaussian frequency window of width a 
centred around frequency A as shown in FIG.|l]), can be 
engineered by choosing the following spatial profile, 

f(x) = n CT e^^ 2aa ( e ^ a + e +jX - 2 ), (11) 

which corresponds to a gaussian multiplied by a super- 
position of planes waves of opposite momentum A. n a is 
a normalisation constant. The spatial profile is therefore 
transformed into the momentum distribution, 
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FIG. 1: 1 + 1 dimensional example of a frequency distribu- 
tion peaked around ±A given by Eq.( 12 1 for a = 1 and A = 5. 



This frequency distribution peaks around around the desired 
frequency A but has a double peaking due to the two expo- 
nential terms. In the 1 + 1 massless case, the field is expanded 
as an integral over ui > and so the peak in the ui < region 
does not contribute. 



This means that in order to couple our detector to a 
peaked gaussian distribution of modes centered around 
A it is necessary to engineer a field-detector coupling 
strength which not only is peaked around the atom's 
trajectory but also oscillates with position. Sharp fre- 
quency modes f(k) = S 3 (k - A) + S 3 (k + A) are obtained 
when f(x) ~ exp(-iA • x) + exp(+iA ■ x). In the mass- 
less 1 + 1-dimensional case the frequency distribution ob- 
tained from a given spatial profile is defined as a function 
of uj > only. Therefore, given the window profile peaks 
are sufficiently narrow and separated, the second peaking 
corresponding to the uj < region, does not contribute to 
the frequency window in this case. The field to which the 
detector couples given by Eq.Q is therefore, 



')-/ 



a u + e 



*4] (13) 



In the general case, the frequency window is peaked 
around two modes corresponding to negative and positive 
momentum. Also see |21) which, complementary to this 
work, looks in detail at the coupling of Unruh-DeWitt 



detectors to Minkwoski modes. Eq.(12) shows that there 
is a trade-off between the width of the frequency window 
and the spatial profile of the detector. 

We now proceed to calculate the instantaneous transi- 
tion rate of the detector model we have introduced. The 
transition rate is a function of the detector's energy gap 
A [9] given by 



dsRe \t 

o 



-iAs 



s)], (14) 



where W(t,t') := (/>(t)</>(t') is the so-called 
Wightman-function and \ip) denotes the state of the field. 
Note that we have assumed that the detector is turned 



on in the distant past. Expanding the field in terms of 
Minkowski modes we find that the vacuum transition rate 
for a stationary detector is 



F T (A) = -9(-A-m)S(A), 



(15) 



where 



S(A) = VA 2 -m 2 |/(-A)| 2 (16) 
and 9 (a;) is the Heavisde theta function defined as 



: x<0 

1 : x>0 



(17) 



Note that in the above result we have explicitly assumed 
that f(k) = /(|fc|) i.e. the Fourier transform of the spa- 
tial profile f(x) depends on the magnitude of k only. The 
result is explicitly independent of the time parameter r. 
This can be traced back to the fact that an inertial trajec- 
tory is a stationary orbit in flat spacetime. The interest- 
ing result here is that the transition rate of the detector 
is tempered by the square of the frequency distribution 
/. We are also interested in analyzing the response of the 
detector when the field has a single Minkowski excitation, 



|1*> 



J d 3 fc$(fc)aj. |0>, 



(18) 



where we define a delta normalised state to have the 
property f d 3 k\Q(k)\ 2 = 5 3 (0) and a properly normalised 
state to have the property f d 3 fc|<I>(fc)| 2 = 1. This state is 
the generalization of Eq.([lO| to three spatial dimensions. 
The Wightman-function in this case is 



W(t,t') := (Ul^r)^')!!*): 



(19) 



Writing the states and the field in terms of the Minkowski 
modes and normal ordering the associated operators we 
find that 



W(t, t') = J d 3 fc|$(fc)| 2 • (0| </»(t)<Kt') 1°) 
+2Re[4(r)/ $ (r')], 



where 



I»(r)= f d 3 k<S>(k)f(k)^- 



(20) 



(21) 



We notice that there are two terms in the Wightman- 
function. The first one corresponds to the vacuum state 
and the second is the contribution from the particle 
present in the field. 

The single particle contribution factorizes into two in- 
dependent functions of r and t'. This allows us to anal- 
yse the transition rate with relative ease. Substituting 
Eq. pOl) into Eq. (fTil we obtain 



tr(A) := f 
Jo 



dse^Uij-s). 



(22) 
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This integral is a Fourier transform of I$(t - s) in the 
s variable and can be computed, either analytically or 
numerically, for specifically chosen / and g. Employing 
the Riemann-Lebesgue lemma, which can only be used 
for functions which are integrable on the real line, one 
shows that /$(t) -»■ as t -> ±00 as long as / and w 
are well-behaved. /$(t) vanishes in the distant past and 
future where the detector is responding only to vacuum 
fluctuations. In other words, the detector only observes 
a constant spectrum in these asymptomatic regions. In 
the intermediate regions the oscillatory response is due 
to the presence of the particle. 



III. ACCELERATED TRAJECTORY 

We now consider a detector following an uniformly ac- 
celerated trajectory. Conformally flat Rinder coordinates 
£ = ( r i £) = (Pi ?i Vi z ) are a convenient choice in this case. 
The transformation between Rindler and Minkowski co- 
ordinates is given by 



and commutation relations satisfy 



t 
x 

y 



a 1 e a ^sinh(ap) 
a _1 e a * cosh(ap) 

y 

z 



(23) 



where a is a postive parameter. This transformation 
holds for the spacetime region \t\ > x which is called 
right Rindler wedge. The coordinate transformation for 
\t\ > -x (left region) differs from Eq. (23) only by an 
overall sign in the x coordinate. The coordinates are 
tailored specifically to the trajectory £ = so that an 
observer travelling along this worldline will measure a 
proper acceleration yJ-A^A^ = a and the proper time 
is parametrised by the coordinate time p. The Klein- 
Gordon equation for a massive bosonic field in a (3+1)- 
dimensional flat spacetime in this case takes the form 

d pp <\> - [% + e 2a t(d yy + d zz )- toV q «] = 0, (24) 
and the solutions are the Rindler modes [HI [TH] 

%,£ 1>a (pi) == N n/a K m/a (Ma-V«)e- iA «(^ 



A a(p,l) := apQ -k ± -x ± , 



(25) 



with M = \fk L ■ k L + m? and Nqu is the mode normal- 
isation constant. The functions i*Qn/ a (-R) are modified 
Bessel functions of the second kind. Here x ± := (y, z) and 
k\ := (fcy,fc z ) are position and momentum vectors per- 
pendicular to the direction of acceleration. f2 is strictly 
positive and denotes the dimensionless Rindler frequency 
and a = +1(-1) corresponds to right (left) Rindler re- 
gions, respectively. The canonical orthonormality rela- 
tions for the 3+1 massive field are, 

K^^j^) = W ~ n')S 2 (k ± - k[)6 aa ,, (26) 



i ] = s(n-n')s 2 (k 1 -k[)5 ac 



(27) 



From our coordinate definitions Eq. (23), we choose 
the detector to be travelling in the right Rindler wedge. 
This then implies that our comoving coordinates can be 
parametrised as r = p and Q = £. The field expansion in 
terms of the parametrised Rindler modes is 

4>(t,£) = f dnd%[u n:ki<+ (p,Oa nrkit+ + h.c.] (28) 



Note the left Rindler modes do not appear in Eq[28] as 
the detector is assumed to be moving in the right Rindler 
wedge. The explicit form of the accelerated detectors 
frequency distribution is 



/(ft, k ± ) = J d 3 ie 2a Zf(OK lQ/a (Ma^e^) e 



+ 2 \z j_ ■ £C j_ 



(29) 



The most significant difference between the inertial and 
the accelerated frequency distributions is the appearance 
of a non-trivial metric factor and the Bessel function. 
Note also that for both massless and massive fields, the 
Rindler modes are defined as an integral over O € K + , 
unlike the Minkwoski mode case. Eq.(29) is a Fourier 
transform in the y and z coordinates, however, it is a 
non-standard integral transformation in the £ coordinate. 
Reminiscent of a Hankel transformation, we expect our 
desired properties of arbitrary mode peaking to still hold. 
Using the integral representation of the Modified Bessel 
function for the second kind 



dt 



(sini(i)) 

gi? cosh(i) 



(30) 



valid for fl/a > and R > 0, we can write the frequency 
distribution as a Fourier type integral that takes the form 



f(k) = f dPipffie** 



(31) 



where now k = (il + 5, k yj k z ) and 

iQf a 



( 1M r*'" - r~ (smh(r)) 2in ' a 



(32) 



S is a phase that is acquired from the integral representa- 
tion of the modified Bessel function. This shows that, in 
principal, the standard properties of the Fourier transfor- 
mation can be used to design a detector profile such that 
we obtain a peaked distribution in momentum space. For 
a concrete example, we shall consider the massless 1 + 1 
field case. The appropriate transformation, in terms of 
Rindler modes, is given by 



/(fi) = f d£e 2a tf(0< 



(33) 
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and the spatial profile we propose in this case is 



f(0 = N(a)e 



•) (34) 



which includes the conformal metric factor that arises 
from the Rindler coordinate transformation. Here N(cr) 
is a normalization constant. This profile reduces the in- 
tegral transformation ( |33[ ) to a standard Fourier trans- 
formation and the resulting frequency distribution is 



(35) 



Substituting this frequency distribution into Eq.Q, we 
obtain 

0(r) = f dnN n/a e-^~ x -V 2 [e- inT a Q , I + h.c] (36) 

In the limiting case where the acceleration is zero Q ->• lj 
and the spatial profile reduces to the Minkowski profile 



given by Eq.( 11 1 



We would now like to expand the field in terms of Un- 
ruh modes which pay an important role in the literature. 
These modes are given by [TH [TS] , 



cos 



: X ,L := cosh ( r n/a)u n k^_+smh(r n/a ) 



,(37) 



n.k, 



where tanh(ro/ a ) = e _7rn '°. Upon parametrising the 
modes with our accelerated comoving coordinates, i.e. 
{piQ = ( T ,£)> an d noting that the left Rindler modes 
have no support in the right Rindler wedge, we find that 
the Unruh modes reduce to 



u 



Q,k ± ,R 



(r,0 



3h(r 0/a ) 



(r,0 



(38) 



Expanding the field in terms of Unruh modes on the 
right wedge results in Eq. ( 39 ) and the corresponding fre- 



quency window operator is given by Eq.(41) 



(r,D = f <md 2 k L N nla K inla ( — e a ^ [co S h(r n/a ) e -^ k ^ a n k R + sirJi(r Wo )e +in ^^a n Al + h.c] (39) 
0(r) = f dnd 2 k ± N n/ J(nX)[cosh(r n/a )e- mT a n ]li R + S i^ (40) 



In case the field is massless and in 1 + 1-dimcnsions, the 
detector interacts with two peaked distributions corre- 
sponding to right and left Unruh modes respectively, 



') = <Mt) + </>l(t) 



(41) 



where 



i B (r) = / dnN n/a e~i° 2 (~ x ~V 2 [cosh(rn/a)e- inT a QtR + h.c.] 
h(r) = f dnN n/a e-^ 2 ^ 2 [sinh(r 0/a ) e + ^a 0iL + h.c] 

I 



(42) 
(43) 



Note that the effective interaction strength is now 
modulated by hyperbolic trigonometric functions. We 
are currently using this detector model to analyse entan- 
glement extraction of sharp frequency Unruh states. 

As uniform acceleration is also a stationary orbit of 
flat spacetime, we expect a time independent vacuum 
transition rate. Using the parametrised Unruh modes, 
we can calculate the transition rate of the accelerated 
detector. We find that for the field in its vacuum state, 



1 



g2^-A/a _ ^ 



S(A) : 



(44) 



where 

3(A) := f d 2 fc 1 [7V?|/(A)| 2 e(A) (45) 

-iV 2 A |/(-A)| 2 e(-A)] (46) 

with ±A := (±A,fcj,,fc z ) and denotes the appropri- 
ate normalisation for the Rindler modes. We can see 
immediately that the transition rate of the detector is 
the expected thermal distribution, where the tempera- 
ture is inversely proportional to the acceleration param- 
eter a, but again modified by the frequency of the field 
operator. We also note that Eq. (44) satisfies the Kubo- 
Martin-Schwinger condition [351 123] 



G 



*V(A) 



F T (-A). 



(47) 



The transition rate is, as expected, independent of time 
due to the stationarity of the trajectory and the invari- 
ance of the vacuum state. Now we shall analyse the re- 
sponse of our accelerated detector model when the field 
contains a single Unruh particle. In the literature, well 
analysed states of the field correspond to maximally en- 
tangled Bell states, see for example, (TUl HH US] . These 
states contain both the vacuum and a single Unruh par- 
ticle. Our starting point will again be the Wightman- 
function which, for the one particle state, takes the form 



W(t,t') := (l p |0(r)^(r')|l P ) 



(48) 



where is a one Unruh particle state defined as 
|l p ):= / df2d 2 fc ± $(f2, k ± )aJ - |0). Continuing in the ex- 
act same fashion as the Minkowski one particle state we 
find 



W(t,t')= J dM 2 fc 1 |$(fc)| 2 -(O|0(r)0(r')|O) 

+2Re[/;(r)/ p (r , )] ) 



(49) 



where I p (r) := / rfOcZ 2 fc x /(Q,fc i )$(fi,fc i )J7n, J ,(r) with 

Un, P (r) = N nla | sinh(r ^ a)e+ ^. p = L (50) 

An Unruh particle has with it an associated wavefunc- 
tion, defined as f (ff2d 2 A; ± $(fi, k ± )Un . p (i, x), which must 
follow an accelerated trajectory peaked around some 
specified centre in position space. These are, by defi- 
nition, accelerated plane waves. Thus, the accelerated 
detector will probe the accelerated plane waves as they 
approach. As the particles wavefunction oscillates as a 
function of r, the detector will observe different phases 
at different times. It is these oscillations that contribute 
to the undulatory behaviour of the detectors transition 
rate. For the Unruh state, the corresponding acceler- 
ated expression for Eq. 22 again has a time dependent 



as for the Minkwoski particle. The Riemann-Lebesgue 
lemma can be used to show that in the asymptotic past 
and future, the response of the detector is the same as 
the vacuum and hence has a thermal signature. As with 
the Minkowski particle, the intermediate regions between 
past and future asymptotic times give rise to an oscilla- 
tory response function. In the limit of high acceleration, 
the second terms in Eq. (49) become negligible and the 



oscillatory integral. It is clear for appropriately behaved 
functions / and g the same analysis can be applied here 



state tends to the maximally mixed state. Thus, the sin- 
gle particle state of the field is dominated by the vacuum 
fluctuations. 



IV. CONCLUSIONS 

We introduce a detector model which naturally 
couples to peaked frequency distributions of Minkowski, 
Unruh and Rindler modes. This detector model is 
suitable for studies of entanglement extraction in non- 
inertial frames. In the (3 + 1) -dimensional case, the 
frequency window of the detector peaks around positive 
and negative momentum inducing a double peaking. 
In the (l+l)-dimensional case, frequency distributions 
naturally peak around a single frequency. We obtain 
analytical results for the instantaneous transition rates 
of the detectors undergoing inertial and uniformly 
accelerated motion. In particular, the transition rate 
of the accelerated detector is the expected thermal 
distribution modified by a smearing function that arises 
from the detectors spatial profile. We also showed that 
the well studied single Unruh particle states produce an 
oscillatory response that is only thermal in the asymp- 
totic past and future. Since for accelerated detectors 
thermal noise is observed in both the vacuum and the 
single Unruh particle state, entanglement is expected 
to be degraded by the Unruh effect for global Unruh 
modes. We also see that the degradation effects occur 
for properly normalised wave packets. 
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